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I. INTRODUCTION 



This article is part of a research program the aim of which is to identify all third-order 
superintegrable systems in two-dimensional Euclidean space. We recall that a superinte- 
grable system is one that has more integrals of motion than degrees of freedom. We consider 
a classical or quantum Hamiltonian 

H=^(p 2 1 +p 2 2 ) + V(x 1 ,x 2 ) (1) 

with two integrals of motion 

X a = fa,jk(x 1 , X 2 )p{p2, a =1,2, (2) 

0<j+k<n 

where pi,P2 are components of the momentum p. In classical mechanics, the integrals Xi i2 
Poisson commute with H, are well defined functions on phase space and the three functions 
{H, Xi, X 2 } are functionally independent. The functions X\ and X 2 do not Poisson commute 
with each other; instead they generate a non-Abelian algebra, usually a polynomial one. In 
quantum mechanics, H and X a are Hermitian operators in the enveloping algebra of the 
Heisenberg algebra (or some generalization of the enveloping algebra). The operators X a 
Lie commute with H , but not with each other. Instead of functional independence, we 
assume that H , X\ and X 2 are algebraically independent. More specifically, we assume that 
no Jordan polynomial in the operators H, X\ and X 2 vanishes. As indicated in (j2J), we 
assume that X a are polynomials in the momentum. The "order of superintegrability" is the 
highest order of these polynomials. 

The best known superintegrable systems are the Kepler- Coulomb system^ with potential 
V = a/r and the harmonic oscillator—^ with V = ar 2 . As a matter of fact, these are the 
only rotationally invariant superintegrable systems in n-dimensions (n > 2) (in agreement 
with Bertrand's theorem^ ). Both of these systems are quadratically superintegrable in 
that the Laplace- Runge-Lenz vector for V = a/r and the Fradkin (or quadropole) tensor 
for V = ar 2 are both second-order in the momenta. 

Until recently, most studies of superintegrability concentrated on the second-order 
case^^ 1 ^ 1 ^. At least in two-dimensional spaces, second-order superintegrable systems 
are well understood. All such systems in Euclidean spaces, spaces of constant curvature, 
and spaces of non-constant curvature with at least two second-order Killing tensors (Darboux 
spaces) have been classified^— 
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Recently, infinite families of superintegrable systems with integrals of motion of arbitrary 
order have been discovered and investigated^^— ^i2£i22i2£i22iMi2si. 

A systematic search for third-order superintegrable systems was started in 2002, both in 
classical and quantum mechanics^. However, the first (to our knowledge), article on third- 
order integrals of motion was written considerably earlier by DracbA He considered the case 
of one third-order integral of motion (in addition to the Hamiltonian) in flat two-dimensional 
complex space in classical mechanics. He found 10 different complex potentials which allow 
a third-order integral. Later it was shown that 7 of them are actually quadratically super- 
integrable and the third-order integral is reducible, i.e. is the Poisson commutator of two 
second-order integrals^ 1 ^ 



The articles in Ref.'s 
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and 
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were devoted to superintegrable systems with one 
third-order and one first or second-order integral. A first-order integral in E<i (R) exists only 
if the potential is translationally or rotationally invariant, i.e. V = V(x) or V = V(r). 
In the classical case, all such potentials are second-order superintegrable and hence known. 
In quantum mechanics, one class of new superintegrable potentials is obtained^ and is 
expressed in terms of elliptical functions, e.g. 

V = h 2 co 2 sn 2 (ujx,k), (3) 

where u and k are constants and sn(ux, k) is a Jacobian elliptic function^. The existence 
of a second-order integral implies that V(xi,X2) allows separation of variables in Cartesian, 
polar, parabolic or elliptic coordinates. Cartesian and polar coordinates were considered 
earlier— 1 ^. The study provided a number of new superintegrable potentials in the classical 
and, much more interestingly, in the quantum case. Indeed, quantum integrable systems 
with higher-order integrals of motion can be quite different from classical ones; a fact first 
noticed by Hietarintar^ 1 ^. The case of quantum superintegrable systems is much richer 
than that of classical ones. Third-order superintegrability with separation of variables in 
the Schrodinger equation in Cartesian or polar coordinates lead to potentials expressed in 
terms of Painleve transcendents (Pi, Pu and Pjv for Cartesian coordinates, Pyi for polar 
ones). The potentials separable in Cartesian coordinates have been intensively studied in 
both classical and quantum mechanics^"—. 

The purpose of this article is to find all superintegrable systems that allow (at least) 
one third-order integral and a second-order integral that leads to separation in parabolic 



coordinates. 



II. THE DETERMINING EQUATIONS 

Let us now assume that the Hamiltonian allows separation of variables in parabolic co- 
ordinates 

x x = ^(£ 2 - T] 2 ), x 2 = £q. 

The quantum mechanical Hamiltonian has the form 

h 2 ( d 2 d 2 \ 



£2 _|_ ^2 yg^ g-q2 



and there exists a second-order integral of the form 



Y = L 3 p 2 +p 2 L 3 + ^J-j {eWM - 77 2 Wk(0) (6) 



with 



'h ^ ( ^ + £ ^ 

^ 2 1 dx 2 £ 2 + rj 2 \ <9£ dr) 

■ 7-/ x ih f ^ d d 

L 3 = -ih{x 2Pl - x lP2 ) = — If— -Vgjij- 

A third-order integral will then have the form 

X= A m{ l I>PipI} + {9i(xi,x 2 ),Pi} + {g 2 (x 1 ,x 2 ),p 2 }. (7) 

j + k+e=3 

The brackets { , } denote anti-commutators, Ajki are real constants and the functions V, 
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fin Cartesian coordi- 



9i 92 obey the four partial differential equations presented in Ref. 
nates). Here we need the equations in parabolic coordinates. To rewrite them in parabolic 
coordinates, it is convenient to replace the unknown functions gi(xi,x 2 ) and g 2 (x±, x 2 ) by 
Gi(£,ri), G 2 (€,rj), putting 

a (x x)~ -^'V) a(x x) - VG1&V) + &2&V) 
9i{xi,x 2 ) — — — — - , g 2 {x 1 ,x 2 ) — — — — - . 



The four determining equations for the integral (JTj), i.e. the commutativity condition 
[H,X] = 0, can be written as 

Gu-^^ = Me,i), (8) 

G 1>v + G u -2 vG ^f 2 ^h 2 (C,v), (9) 

g 2 ,, + iG t \ ~ ^ 2 = h (£,v), (10) 



GiV^ + G 2 V V = ^-0, (11) 

with 



2F 2 ^+2F 3 K, 

£2 +?? 2 

F 3 yg+3F 4 K, 

£2 +r? 2 ) 



and 



•inn 



,2^2 



(12) 



| (3^ + 2r,F 2 - £F 3 ) % + ( V F 2 - £F 2 - 7]F 3 + £F 4 ) + ( V F 2 - 2£F 3 - 3 V F A ) V m 

- p f -^y f| - fl ' -!w + -fl + + „) vj 

The subscripts denote partial derivatives and the expressions F\, F 2 , F 3 , F4 are polynomials 
in £ and 77: 

Fx = ?? 3 ^oo3 + £v 2 Aoi2 + fA 030 - - rf (f + rj 2 ) A W2 - ^7] 2 (£ 2 + r/ 2 ) A m - - r]f {f + ^120 
+\v 3 (C 2 + r? 2 ) 2 A 201 + \^ 2 (e + n 2 ) 2 A 2lQ - l -^(e + ri 2 f A m » + ev A021 

F 2 = 3 i i 1 2 A 003 - 7] (if - 2 f) 4)12 - 3 frj A 030 - ^ £ V 2 (£ 2 + r/ 2 ) ^102 + ^ 3 (£ 2 + r] 2 ) A U1 
+\ e + 2V 2 ) (f + V 2 ) A l20 - \tr? if + V 2 ) 2 A 201 - \rj ( V 2 + 2£ 2 ) (£ 2 + r, 2 ) 2 A 2W 

+hv 2 (e+v 2 ) 3 A 30Q -A Q21 z (2 V 2 -e) 



F 3 = 3 C 2 V A 00 3 ~ t (2 r? 2 - £ 2 ) A 012 + 3 £ r? 2 A 3o + ^ £ 2 (£ 2 + r/ 2 ) A 102 + {f + r, 2 ) A 1U 

- 1 -v(v 2 +2f)(e+ rf) a 120 -\^e{e+ v 2 ) 2 a 201 +\ae+^ 2 ){e+ v 2 ) 2 a 210 

-l^e (e 2 + r? 2 ) 3 ^300 + ^02i (r? 2 -2^ 2 ) 
o 

F 4 = fA 003 ~ i 2 r] A012 - rfA mo + {f + rj 2 ) A 102 - - r) £ 2 {f + r? 2 ) A m + - £ r? 2 (£ 2 + r/ 2 ) Ai 20 
+| if + V 2 ) 2 A 20 i - {e + r, 2 ) 2 A 210 + ~ f (£ 2 + r? 2 ) 3 A 300 + ^^021. 



The determining equations flBlfTOj) are the same in classical and quantum mechanics but 
ffTTj) contains the Planck constant on the right hand side. The corresponding equation in 
classical mechanics is obtained by taking the limit h — > so (|TT|) is greatly simplified. 
Hence the difference between classical and quantum integrability (and superintegrability) 
for third-order integrals of motion. 

The system flSl fTTj) is overdetermined. The first three equations imply a linear compati- 
bility condition for the potential 

= F 3 V^ + (3 F 4 - 2 F 2 )V^ V + {3F 1 -2 F 3 )V im + F 2 V vm 

+ (2(F 3C - F 2v ) - 3 ^~ff + y^ ) % + ^ _ Fh) _ 3^-s^3 + 7^ Vm 

+ (2 (3F Vn - F 2 , - F 3v + 3F H ) - ^f^P^ ) V, v 

+AVr, + BV(, 

where 

■A - t 2ryr] - At Mi + 3^ 4€f H — + 

+2 2U 2 F 4 + F 2 ^ 2 + 7 r/ 2 F 2 - 12 £ r; F 3 + 3 F 4 r/ 2 



(13) 



-S — 3Fi vv — 2F 2rj £ + 
+2 



(e 2 + ^ 2 ) z 

21r7 F lv + 3£ F u - 6£ F 2?? - 6r/ F 25 + 77 F 3?? + 7£ F 3? 



F 3 r/ 2 + 3 F^ 2 + 21 r) 2 F x - 12 £ 77 F 2 + 7 £ 2 F ; 



£2 + rj 2 
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(e+r/ 2 ) 2 

Compatibility between the first three determining equations (IH tiTOj) and the fourth one 
ffTTj) requires three more conditions, this time nonlinear ones. Indeed, solving ffTTj) for G 2 , 
we have 

G2 = y (j® - , V„?0. (14) 



Replacing G 2 from (1T41 into (IHl fTOj) . the system can then be solved for G±, 




G 2 



h 5 h 2 ^ - AV^h 4 



(15) 



with 



h, = (e + rj 2 ) (Ah 3 v v 3 + Ah 2 v^v 2 + - h 2 v 2 % - h 2 v^v v ^ - h 2 (v m v v - v^) $) 

+A V V 2 - AryV 2 - 8£V v Vt 
h = 4(e 2 + v 2 ) (VMiYa - v m + (V 2 - V 2 )V^ + 4faV € - £V V )(V 2 + V 2 ). 



Replacing (fl~5l) . into (IHl fTOj) . gives the three additional non-linear compatibility conditions on 
the potential, namely 



\hj v V MM, ) i 2h 5 V v (e + V 2 ) 

V ; „ ^(e + r/ 2 ) 

III. GENERAL FORMS FOR Wi AND W 2 



In order to determine all possible potentials which separate in parabolic coordinates and 
admit a third-order integral of motion, we begin with the linear compatibility condition 
(JT3J) • Replacing V with the form as in §5§, the compatibility condition can be differentiated 
to obtain a system of linear ordinary differential equations (ODEs) for W\ and W 2 . These 
equations are given in Appendix |A] for W\. Notice that interchanging £ and rj has the effect 
of changing the sign of the coefficients to (—iy +k Ajke- Thus, equations for W 2 are of the 
same form as (1A1HA141) . up to a change in sign in these constants. We begin our search for 
the admissible potentials by solving equations (IAlHA14p . 

Theorem 1 Given a Hamiltonian which admits a third-order integral of motion with a 
potential which separates in parabolic coordinates as in |3p. Then any admissible terms in 
the potential are included in 




(16) 



16 



o^ln (f + y/e + #i) + a 2 £ 



(19) 



j=0 
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(20) 



/3t In (rj + ^fTB^ + /3 2 r] 

In particular, any solution of the system ( Fl3l) has the form (jSJ) with W\ and W2 as in ([19 
and (120]) with appropriately chosen constants. 

Proof: Beginning with (lAlj) . in the case that v4 30 o 7^ 0, the solutions are given by 

14 



^i = 5^Qf + c _ 2 r 2 + c- 4 r 4 - 



(21) 



i=0 



Similarly, for flA2j) . in the case that A 2W 7^ 0, the solutions are 

15 



^i = $>e + c_ 2 r 



(22) 



i=0 



Equation (1A3I) in addition requires that C14 = ci 5 = 0. 

Now assume A 300 = A 210 = 0. This sets (IAHIA3I) to be satisfied identically and fl A4[) 
becomes 



/ d 17 d w d 15 \ 

- ( (f ^201 + 2A 120 ) — + 33£ A 201 — + 255A 201 — J W 1 = 0. 

The solutions of ( 1231) . assuming A 20 i 7^ 0, are 

w ai f l n + v^T^T) + a 2 £ 

with A\2o = B1A201/2. If A 20 i = 0, the solutions are 

16 



(23) 



(24) 



(25) 



i=0 



Now assume A 300 = A 2 io = A 201 = Ai 20 = 0, ()Al-HA4j) are then identically satisfied and ( 1A5j) 
becomes 



/ d 17 d 15 d w \ 

i(A W2 e + 2A 02 i) -—j + 255A 102 — + 33£A 102 — J W x = 0. 



(26) 



de 7 de b dt 16 

The solutions of ( 1261) for Ai 02 7^ are as (l24"j) with A 02 i = 2A W2 Bi and for Ai 02 = are as 

(125}. 

If we now assume A 300 = A 2W = A 201 = A 120 = A W2 = A 021 = 0, (fA6l) becomes 

d 17 d 16 d 15 ~ 



A 



IKK! 



3^3717 + 111^ + 969 



de 6 di 



15 



Wi = 0, 



(27) 
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the solutions of which are given by (12T1) . 

Now, assume A 3Q0 = A 2W = A 20 i = A 120 = A 102 = A 021 = A 003 = 0, this implies (TM]- 
IA6[) as well as (IA10HA121 are identically satisfied. Equation flA7h becomes 

UAme + QA 030 -AA 012 ) <F„, / 51 W 16 



(/(l7 ^i+ ( Y^iur- 36 A 12 + 54A 030 J 

813f Am d 15 „ d 14 

Wl + 2016A m -— = 0. (28) 



2 d£ 15 1 111 d£ 14 

The general solutions of (|28|) are of the form 
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W 1 = £ + c_ 2 r 2 + V J^== 1 , (29) 

with A 030 = 1-BiAui + |v4 i2 for A m ^ and of the form ([22]) for Am = 0. 

Now assume that A 300 = A 2W = A 201 = A 120 = A W2 = A 021 = A 003 = A U1 = 3A 030 - 
2A i 2 = 0, this implies that (IA14j) has the form 



/ d 17 d 16 d 15 d u \ 

A u + 57£ 2 ^ + 1023^ + 5760^J W 1 = 0. (30) 

The solutions of (I3"0"j) for A 012 ^ are given by 

14 

w, = Y, a? + c- 2 r 2 + c_ 4 r 4 + c_ 6 r 6 - (31) 

i=0 

This is the final case, since assuming A i 2 = in addition to A 300 = A 2W = A 20 i = Ai 2Q = 
A\o 2 = A 021 = A 003 = Am = 3A 030 — 2A i 2 = has all the A's equal and so there is no 
longer a third order term in our constant of motion. 

Thus, the most general form of the function W\ is given by (fT9|) . By direct analogy, the 
most general form of the function W 2 is given by (T20l) . 

□ 

Corollary 1 The potential for any 3rd-order superintegrable system which separates in 
parabolic coordinates satisfies a non-trivial system of linear ODEs for both W\ and W 2 . 

Proof As shown above, the compatibility conditions ( fl3|) are satisfied identically if and 
only if all the A^ are 0. On the other hand, if there is a non-zero A^k then the functions 
W\ and W 2 in the potential will satisfy some linear ODEs. 

□ 
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IV. THE ABSENCE OF IRRATIONAL TERMS 



In this section, we show that the only possible form of the potential IS cLS cL rational 
function of £ and rj. We shall show this by contradiction. Namely, we consider the case that 
either a>i 7^ or both a 2 7^ and B x 7^ in (|TU|) . 

If we substitute the general form ( |T9|) into equations (1A1HA14]) . it is immediate from the 
previous section that A^ Q0 = A 2W = 0. We also obtain the following possible restrictions on 
the constants, as suggested in the previous section: 

Bi(A 10 2-A 12 o) 



A-003 — A021 

A012 = A, 

^201 



2 

_ Amgi 
030 — 2 ' 

2A120 



B 



(32) 



In the case that B\ ^ 0, W\ becomes 

Wx = c + cxC + c 2 £ 2 + c 3 £ 3 + c 4 £ 4 + c 5 £ 5 + c 6 £ 6 + 
with the following cases 



C-2 



+ 



7F+sT 



, (33) 



{Bi ^ 0, c_ 2 = C5 = c 6 = 0} 
{Si ^ 0, A 201 = 0, c_ 2 = c 6 = 0} 

{Si ^o,A m = ^201 = o,c_ 2 = 0} 

{5x^0,^02 = ^201=0.} 

In the case that a\ 7^ and Sx = 0, Wx becomes 

Wx = c + C i£ + C2 e + C3 e 3 + c 4 e 4 + c,e + c Q e + ^ + 

with the following cases 

{«i 7^ 0, Si = 0, c„ 4 = c 5 = c 6 = 0} 
{ai ^ 0, Si = 0, A201 = 0, c_ 4 = c 6 = 0} 

0} 



(34) 



+ a<i ln(£), 



(35) 



(36) 



{ai 7^ 0, Si = 0, A m = A201 = 0, c_4 

{«1^0,Si = 0,^102=^01 = 0}. 

Note that when Si = 0, a 2 becomes an additive constant which is absorbed into Co. 

Next, we substitute these cases f[3"4"|) or ([35]) . along with the forms of W\ fl3"31 or (1331) and 
W2 (POL into the compatibility condition ( fT3l) to obtain 

(37) 



= M 2 /n (V + V^ 2 + S 2 ) + M lV V + S 2 + M 
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where the M; are rational functions which are too long to be presented here but are available 
from the authors upon request. Solving M 2 = gives 3 possibilities: either (3\ = 0, B 2 = —B\ 
or, in the case that B\ is not zero, B 2 may vanish if additionally 4 m = and A 102 = 
BiA 201 /2. 

Solving Mi = and M = gives similar forms for W 2 '■ If B\ 7^ 0, then W 2 has the form 

-air] In ( n + ^n 2 + B 2 ) + (3 2 n 
W 2 = d lV + c 2 n 2 - C477 4 + c 6 r/ 6 + -± + ^ — — ^ (38) 

and if -Bi = 0, then it has the form 

W 2 = c 2 r, 2 - c 4 7/ 4 + c eV 6 + ^l-t±- ai Info). (39) 

In both cases, C3 and C5 are also required to be as well and as several additional sets of 
constraints which are required to completely solve (fl3l . We shall return to some of these 
cases later. 

To obtain a contradiction for the potentials admitting logarithmic singularities, we now 
turn our attention to the non-linear compatibility conditions. Beginning with Eq. ( Tl6|) . we 
clear the denominator and consider the equation 

'hA r]h 2 h 5 ^ + 4h 4 (£V v + riVt) 



o = hiv v (e+v 2 



(40) 



hsJz Ah 5 V v (e + V 2 ) 

In the case that both a\ and B\ are assumed non-zero, equation ( 140]) will have polynomial 
dependence on the quantity In (^ 2 + a/£ 2 + B\ j . Substituting the forms of the potential 
obtained above ( 133]) and ( 138]) into this quantity. ( 140]) . gives 

= 32a 6 B 3 1 Kln 6 (f 2 + v^ 2 + K) + O (in 5 + v^ + ^i)) , (41) 

^ = 9A m ^ 14 + ^(16^ + 550^ 

+ (2(48A 102 - 23A 201 ) V 2 - B x {l§A im - 1L4 201 )) r^ 11 + 0(£ 10 ). (42) 

Therefore, since it was assumed that «i 7^ and B\ 7^ 0, the condition ( 1401) requires that 
An = A02 = ^201 = which is a contradiction because in this case all of the A^'s would 
be identically zero. In the case that B\ = 0, (14*0]) is a polynomial in ln£ with leading order 
term 

0= (eTWW ln4e + o(ln3 ° (43) 
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K = 192a^ 111 ^ 12 r ? 2 - 640c^ 102 £ i y - 32(864d_ 2 c 4 A 102 + a 2 A 201 )a^ 13 rf + ... (44) 

For simplicity, we give only the most relevant terms instead of the highest order ones in K. 
From these terms, it can be seen that the condition (140 j) leads to a contradiction since they 
would imply A201, ^102 and Am are all and so every Ajki would vanish. Thus, there are 
no logarithmic singularities in the potential. 

Next, we turn in particular to the case that a± = /?i = 0. In this case, we proceed 
slightly differently because computationally, it is more difficult for MAPLE to compute 
coefficients of the expression (140p with respect to (£ 2 + rf) k ^ 2 without first simplifying the 
entire expression. On the other hand, unlike in the previous section, we can solve (151 1TUj) 
without much difficultly and replace the integrated forms of G\ and G 2 into ffTTj) to obtain 
the needed contradictions. 

In this case, we require the complete solutions for (1T5|) . Namely, W\ must satisfy (155|) 
with C3 = C5 = and W 2 must satisfy ( 1551) . both with a% = f3\ = 0. Additionally, the 
constants satisfy one of the following cases: 

{B 2 = -B t , A 20 i = A 102 = 0, c = ci = di = c 6 = 0} (45) 

{B 2 = -B x , A 201 = A U1 = 0, c = ci = di = c_ 2 = 0} (46) 

{B 2 = -B u A m = 0, A 102 = A 120 , Cl = d x = c 6 = c_ 2 = 0} (47) 

{B 2 = 0, Am = 0, A 102 = A m , c = c 1 = d 1 = c 6 = c_ 2 = 0} (48) 

{B 2 = -B u Co = c 1 = d 1 = c 6 = c_ 2 = 0} . (49) 

There is also a complex solution for the constants 

{B 2 = -Si, A 20 i = 0, Am = iA W2 , d x = ic h c = c 6 = c_ 2 = 0} . (50) 

To obtain the needed contradictions, we use the obtained sets of solutions for (fT51) to 
solve (151 1T0]) for G\ and G 2 and replace these solutions into ( TTTT) . For example, in the case 
identified in (1431) the relevant G's are 



c __ , '(-SB.e + ^e^ + ^B^a, £77 (-r] 2 + 2 B x + j 2 ) 

2^/eTBi{e + v 2 ) Vv^b; (e + v 2 ) 

1/2 (2 ^ + V 2 Bl + 4 - 3 B,e) t c 4 - fr 2 "* 2 "^ " 2 ) A m 
^2 + ^ 



fel + ^2 + ^3 (51) 
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c = ^ v (- V 2 + 2B 1 + e) 02 1/2 (-Bie + 4 £ V + 3 V 2 B,) ft 

+1/2 77 (2 - B,e + 4 e V + 3 rfB,) c 4 + ^ ~ ^ A x 



ii 



( . -A* - ^7^3 + & 2 , (52) 

where fci,A;2 and are constants of integration. These solutions for G\ and G 2 when 
substituted into ( ITT1) give 



= T lV /e 2 + 5i vV - #i + T 2 ^e + Bi + Tgv 7 ^ 2 - #i + T 4 . (53) 

Again, the Tj's are rational functions of £ and 77, and can be obtained from the authors. 
Solving these systems, we obtain (3 2 = and a 2 = or Am =0 , which gives a contradiction, 
since in these cases either all of the Ajke are zero or the potential reduces to a rational 
function. By checking each case in this manner, we find that when B\ is assumed to be 
non-zero the potentials reduce to rational functions. 

Thus, we have shown by contradiction that the only possible potentials which satisfy 
both the linear ffl3|) and nonlinear (I16H18P compatibility conditions are rational functions of 
£ and rj. 



V. FINAL LIST OF SUPERINTEGRABLE POTENTIALS 

Since there are no longer any irrational terms in the potential, it is a straightforward 
computation to find the admissible choices of constants which satisfy the linear compatibility 
condition ( 1131) . to use these choices to solve the linear partial differential equations (IBl TIOl) for 
G\ and G 2 and to solve the resulting algebraic system determined by the coefficients of (1TTI) . 
In this section, we exhibit the possible potentials which remain, i.e. those potentials which 
separate in parabolic coordinates and admit a third order integral of motion. Remarkably, 
the only such potentials are second-order superintegrable. These are exactly the potentials 
which separate in parabolic coordinates as well as another orthogonal coordinate system in 
E 2 (M). We also obtain a potential in E(l, 1) which is presented in Appendix iBl 

It is interesting to note that, in addition to the systems which admit a single third- 
order integral, this method allows us to obtain potentials in the quantum case which admit 
more than one third-order integral. Furthermore, these additional potentials give proof 
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that the quantum correction to (TTTj) . namely $, is not identically for potentials which 
are superintegrable in both the classical and quantum cases. The quantity $ only vanishes 
when, in addition, the appropriate choices of A^i are assumed. For example, for potential 
Vi below, the quantity $ vanishes only when all of the A^'s are assumed except for A m2 . 



A. Potentials which admit a single third-order integral 

1. A deformation of the anisotropic oscillator potential: V\ 

The following potential 



a(Ax 2 + y 2 ) + 2fix+ n2 



Vi = (v A ~ £V + e 4 ) a + m 2 ~ V 2 ) + ^ (54) 

7 



ar 



v 

(3 cos 2 9 + 1) +2(3 r cos 9 H — 9— 

v sin 9 



admits a third-order integral with all all Ajke = except A021 and functions 

a 012 (-2 « ev + v 4 a e 4 + ^ e 2 + 7) 



d = 2 
G 2 = -2 



^2 

A12 (-2 r/ 2 a e 6 + £ 4 - ?/ 2 /3 £ 4 + 7) 



77 £ 

2. A deformation of the Coulomb potential: V2 

The following potential 

t 2 + V 2 V£ 2 ' V 



a /3 

+ 7 



2\J x 2 + y 2 \ x + a/x 2 + y 2 x — y/ 'x 2 + y 2 

a 8 7 



2r 2 (cos + 1) 2r 2 (cos 0-1) 2r 
admits one third-order constant of motion with all Ajke — except A210 and functions 

A 2W (2 v A l e 2 + 2 r/ 4 a + 4 £ 2 /3 ?? 2 - n \ 2 h 2 + 2 £ 4 /3) 



Gi 
G 2 = 



4^ 2 

A10 (2 £ 4 /? + 2 £ 4 7 ?? 2 + 4 £ 2 a r? 2 - r? 2 ^ 2 + 2 ry 4 q) 

477 e 2 
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3. A second deformation of the Coulomb potential: V 3 



The following potential 



V * = C2 i 2 56 

£2 _|_ ^ 



= » / ! (oi\/x+ \fx 2 + y 2 + (3 J ^x 2 + y 2 -x + 7) 

1 ( 6 6 \ 

= — I av^cos - + f3y/2 sin - + 7 I 
2r \ 2 2 J 

admits one third-order integral associated with A 102 (all the remaining A^t = 0) 

Gi = ~A 102 (-/3^ + /3»7 2 + 2J7ae + 2J77) 
G 2 = ^A W2 ( a ^ 2 -ar ] 2 + 2^ 1 + 2r ] ^/3). 



B. Potentials which admit more than one third-order integrals 

The following potentials are sub-cases of those from the previous section, in the quantum 
cases. They admit at least two third-order integrals. 

1. V 1 subcases 

The following potential 



v ha = ne- v 2 ) + it 2 (57) 



h 2 

h 2 

2/3* + %. 



y 2 

admits three linearly independent third-order integrals associated with the constants 
A i 2 ,A 0Q3 and A W2 (the remaining Ajkgs are 0). The functions G\ and G 2 are 

(3h 2 V 2 -h 2 e + 2y 4 ef3)A l02 o h 2 A 003 n ( V *ef3 + h 2 )A 0l2 
(-3h 2 e + h 2 V 2 + 2ev^)A 102 o h 2 A 003 n (ev 2 P-h 2 )A 012 

L ' 2 ~ ^71 1" 6 c 2 V 1 C2 " 

2^2 ^2 V £2 
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2. V 2 subcases 



The following potential 



£2 + ^2 ^2^2 

7 /I 2 
= tt- + 



2r r 2 sin 2 9 ' 

admits two linearly independent third-order integrals associated with the constants A 300 and 
A210 (the remaining Ajm's are 0). The functions G\ and G 2 are 

fr 2 (£ 2 + ?? 2 ) (3 e 4 + 2 gV + 3 r? 4 ) A 3 qo (3 r? 2 ^ 2 + 2 t? 4 7 , £ 2 + 2 rftf + 2 £ 4 fr 2 ) A 210 

/i 2 (e 2 + y 2 ) (3 e 4 + 2 ey + 3 t? 4 ) a 3 qo (3 77 2 e 2 /^ 2 + 2 e 4 /* 2 + 2 t^ 4 7 + 2 r? 4 ^) a 210 

8£ r? 2 4?7 £ 2 

Another subcase of the potential V 2 which admits three linearly independent third-order 

integrals is given by 

7 ft 2 



G 2 = 



2r 2r 2 (l + cos#)' 

The integrals are associated with coefficients A 30 o, ^210 an d ^201 with the remaining Aj k £ = 
and functions 

^(3r ? 2 + 2a(e 2 + ^ 2 )^300 , (2 /^V + 2r ? 2 e 2 7-n 2 )^io 
^ (2e 4 7-3/iV) 



G 2 



4^2 

/i 2 (3 ?? 2 + 2 £ 2 ) (£ 2 + V 2 ) A 300 ?7 (2 e 4 7 + 2 /*V + 3 h 2 e) A 



210 



8£ 4f 
A201 (2£ 4 7 -/*V) 



4£ 

C. Potentials admitting group symmetry 

The following potentials admit a Killing vector associated with group symmetry. 

1. Rotational symmetry 

The Coulomb potential 

ot ct , . 

£, +V r 
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admits four linearly independent third-order integrals associated with constants -4 30 o, -4 2 io, ^201 
and A W2 as well as a Killing vector associated with fa in 

£(2ar ] 2 -h 2 ) r ] (2ae + ^) v (£ 2 + V 2 ) 
k-i = ^210 ^ V A 102® - A 20l fa 

Cr 2 = -A210 : h -4io2^« + ^201 : h «i ~ • 



2. Translation symmetry 

The first potential is 

^- = w = 7- (61) 

It admits 4 linearly independent constants of the motion associated with the constants 
A 012 , A 030 , Am, A 2 \o and a Killing vector associated with fa, with functions 

n a (rf + e 4 ) ^210 a (ry - Q (ry + Q 4i A^a 
Gl = v 1 c 2 

a (r? 4 + e 4 ) -4 210 a (77 - Q (77 + Q A m A 012 q 

~~ W 2 — + ^ 2 -^~ hri - 

The potential 

^2 £2 

^,v,» = ^ = -2 ( 62 ) 

admits 8 linearly independent third-order constants of the motion associated with constants 
A003, -4oi2, -4 3o, -4io2, Am, A 20 i, A 2W , and A 300 as well as a Killing vector associated with 
the constant fa. The G"s are 

_ fr 2 (3 r? 4 - e 4 ) -4 201 fr 2 + e 4 ) Agio fe 2 (3 rf + 5 £ V + 5 g V + 3 £ 6 ) -4 300 



^ 2 (3e 4 -^ 4 )-4 20 i ^ (V 4 + e 4 ) ^210 /i 2 (3?? 6 + 5£y + 5£V + 3£ 6 )A 



s~1 • v \^ 'I / -* *ZU1 V'/ 1 "S / ^ziu . • v \^ • 1 ' ^ ^ '/ 1 ^ ^ '/ 1 ^ ^ j * ^300 7 

2 ~ i^ 2 ^e 2 w 

Finally, the potential 

Vt,* = «(e 2 - V 2 ) = 2«x (63) 

admits four linearly independent third-order integrals associated with constants -4 00 3, A W2 , 
A 021 , A i 2 and a Killing vector associated with fa, with 

Gi = -r] 3 £ 2 A 102 a - 2 r] 3 A 021 a + 2 77^ A m2 a + 2 77 £ 2 A 02 ia + 77 fa 
G 2 = -£ 3 A 102 a, r] 2 + 2 £ 3 A Q2l a + 2 77 A m2 £ 2 a - 2 £ A 21 a, rj 2 + £ fa. 



17 



VI. CONCLUSIONS 



The main conclusion of this article is a negative one: all third-order integrals for potentials 
separating in parabolic coordinates in E 2 (W) are reducible. The corresponding potentials 
allow separation of variables in at least two coordinate systems and are already known to be 
quadratically superintegrable^" 3 ^. Thus V\ (|54j) is a deformed anisotropic harmonic oscillator, 
separable also in Cartesian coordinates. The potential V 2 fl55l) is a deformed Coulomb 
potential, separable also in polar coordinates. The potential V 3 (|56|) is a different deformation 
of the Coulomb potential, separable in two different parabolic coordinate systems (with the 
directrix of the parabolas either along the x or the y axis). Since these potentials are 
second-order superintegrable, they are the same in classical and quantum mechanics. 

This is in stark contrast with the results obtained in the case of potentials separating in 
Cartesian^ or polar coordinates^. There the results were much richer, mainly because of 
the role played by the linear compatibility condition ([131) . 

Indeed, let us first consider a potential separating in Cartesian coordinates V = W\ (x) + 
W 2 (y)- The linear compatibility conditions for this potential, written in Cartesian coordi- 
nates is given by 

— F 3 Wi, xxx + 2{F 2 ^ — F 3iX )Wi iXX — (3Fi m — 2F 2 , xy + F 3tXX )Wi iX = 

F 2 W 2my + (F 2j2/ - F^ x )W 2m + [F 2m - 2F 3<xy + 3F^ xx )W 2y (64) 

with 

Fx = A 30 o2/ 3 + A 2W y 2 - A 120 y + A 030 

F 2 = 3 A 300 xy 2 - 2 A 2W xy + A 201 y 2 + A 120 x - A ul y + A 021 
F 3 = -3 A 300 x 2 y + A 2W x 2 - 2 A 201 xy + A in x - A 102 y + A 012 
Fa = A 300 x 3 + A 201 x 2 + A 102 x + A 003 . 

Differentiating (I641) twice with respect to x gives two linear ODEs for W\. 

(3 A 3m x 2 + 2 A 201 x + A 102 ) W[ b) + (36 A 300 x + 12 A 2m ) w[ A) + 84 A 300 W® = 0, (65) 
{-A ul x - A 210 x 2 - A 012 ) Wi 5) + (-12 A 2W x - 6 A in ) - 28 A 2W W^ = 0- (66) 

Differentiating (1641) twice with respect to y gives two linear ODEs for W 2 : 

(3 A 3m y 2 - 2 A 2W y + A 120 ) W 2 {5) + (36 A 300 y - 12 A 2W ) W 2 {i) + 90 A 3m w[ z) = (67) 
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(A 201 y 2 - A in y + A 021 ) W 2 (5) + (12 A 201 y - 6 A m ) W 2 (i) + 30 A 201 W 2 {3) = 0. (68) 

Interesting potentials in quantum mechanics, involving elliptic functions or Painleve 
trascendents are obtained if either of these pairs of linear equations are satisfied trivially, i.e. 
if the coefficients in (I65H66P or (I67H68P vanish identically. Both systems of equations (I65H68P 
vanish identically if and only if the only non-zero coefficients Aj^g are A 030 and A 003 . The 
potential in this case is for instance 10 

V(x,y) = h 2 [(J^Pjiuix) +ulP I (u 1 y)] , 

where Pi is the first Painleve transcendent^ 4 -. If, for example, ( 167H681) are satisfied trivially 
but ( 16311661) are not, this leads to potentials of the form 1 ^ 

'46\ 1/3 \ „ / 4b 



V(x,y) = ay + (2h 2 b 2 ) 1/3 



where Pjj is the second Painleve transcendent^ 4 -. 

Now, let us consider potentials allowing separation of variables in polar coordinates 3 ^ 

V = R(r) + ^S{6). (69) 

The linear compatibility condition reduces to 

= r*F 3 R"' + (2r 4 F 3 , r - 2r 2 F 2fi + 3r (2 r 2 F 3 - F x )) R" 
+ (r 4 F 3 , rr + 2 r 2 (3F 3 , r + 3F 3 - F 2 ^ e ) - 4 rF 2>8 + 3F 1<ee ) R' 

+^F 2 S>" - 4, (2r 3 F 3 , r - 2rF 2fi + 6F a ) S" , (70) 

+ 73 (3r 5 F 4 , rr + 6r 4 F 4 , r - 2r 3 F 3 , r0 + 3r 2 F %T + r(F 2)6e - 2F 2 ) - \2F lfi ) S' 
(2r 4 F 3 , rr - 12r 3 F 3 , r + 4r 2 (3F 3 , r - F 2 ,. e ) + 4rF 2fi + QF lfi e + 18Fi) S 
where here the F^s are given by 

F x = Ax cos 36 + A 2 sin 36 + A 3 cos 6 + A A sin 6 

F 2 = ~ 3Al Sin 3e+3Aa COS r 3e -^3 cos 9+A 4 sin 9 + ^ CQs26 + sin2Q + Bq 
p _ -3Ai cos 36>-3A 2 sin 39+ A 3 cos 9+A 4 sin 9 , -2Bi sin 26+2B 2 cos 29 , q ^ Q \ (J g j n Q 
7Ti _ Ai sin 39— A 2 cos 39- A 3 sin 9+A 4 cos 6 Bi cos 29+B 2 sin 20+Bn Ci sin0-C2 cos 6» , n 

-^4 — 13 15 ; r J-Sq, 



with 



Ax 


Ao30— A012 

4 ' 


A 2 


A, 


3v4oo3+A 21 

4 ' 


Bx 


Cx 


= ^210; 


c 2 



A021 — Aqqz A 3Ao3Q+Aqi2 

4 > ^3 — 4 > 

2 ) n 2 — ~2 1 n ° ~ 2 ' 

^201) D = A 300 . 
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Equation ( ITU]) is satisfied trivially if all the Ajki = except A 300 . This leads to the potential^ 6 - 



where R(r) is arbitrary and P(8,t 2 ,t 3 ) is the Weierstrass elliptic function 3 . This potential 
allows a third-order integral, however it is algebraically related to the second-order one and 
the system is hence not superintegrable. If however we consider the subcase R(r) = 0, 
Eq. ( JT01 simplifies. It being satisfied trivially allows further constants in the integral to be 
nonzero, namely A 300 , A 210 and A 2 oi- This leads to a superintegrable potential expressed in 
terms of the sixth Painleve transcendent^, Py/(sin#/2). 

In section IIIH we have shown that the linear compatibility condition (fT3l) is never satisfied 
trivially for a potential separating in parabolic coordinates. This rules out any dependence 
of W\ and W 2 on elliptic functions or Painleve transcendents. 

In addition to the real potentials presented in section |V] we have recovered a complex one 
(see Appendix [Bj which is known to be superintegrable 22 . It can be transformed into a real 
form; however, it will live on the pseudo- Euclidean plane E(l, 1) rather than the Euclidean 
one. 
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Appendix A: System of linear ODEs for W\ and W 2 

The linear compatibility condition f lT3|) . implies (by differentiation) two systems of 14 
ODEs for W\ and W 2 , respectively. Those for W\ are given below and those for W 2 can be 
obtained from these under the interchange £ — > i] and Aj^e — > (—iy +k Ajk£. 



V(r,9) = R(r) + -P(e,h,t 3 ), 




(Al) 
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0= A 2W (tw[ 17) + 18Wf 6) ) (A2) 

= A 2W (ll fwi 17) + 503£ 2 W[ 1&) + 7145 £wf 5) + 31360wf 14) ) (A3) 

= -9^300 (VvFi (17) + 65£ 3 VFi (16) + 1485£ 2 Wi (15) + 14070^i (14) + 46410Wi (13)N ) 

v ' (A4) 

+ (4oi£ 2 - 2A 120 )wi 17) - 33£4oi W[ m) - 2554oi w[ lb) . 

= -f 4oo (^W^ 17 ) + 93£ 5 W^ 16 ) + 3375£ 4 W/ 15 ) + 60900£ 3 VF/ 14 ) + 573300£ 2 W^ 13 ) 
+2653560£Wi (12) +4684680 Wi (11) ) 

_A 201 (^V/ 17 ) + ^ + ^ W 15 ) + 27615^i {14) + 83265 W^) (A5) 

(b^Wx^ + 169^i( 16 ) + nmwi^ 



+ (£ 2 4o2 + 2^021) Wx 17 + 33£4o2 Wi (16) + 2554o 2 W / i {15) . 
= =^4oo (e 8 ^i (17) + 12ie 7 VFi( 16 ) + 5993^ 6 VFi( 15 ) + 157962£ 5 Wi( 14 ) + 2410590£ 4 W^ 13 ) 
+21676200£ 3 VFi (12) + 111351240£ 2 Wi (11) + 296215920^1 (10) + 309188880 Wi (9)X 

_A 201 (5f + ^ Wl (16) + 1483111 ^(15) + 124418 £3 ^(14) 

+1081626£ 2 Wi (13) + 4585308£VFi (12) + 7339332 Wi (11) ) 

-A120 (7£ 4 VFi (17) + 388C 3 VFi (16) + 7464^ 2 Wi 15 + 58632^i (14) + 157248 Wi (13) ) 
+2A 102 (£ 4 Wi (17) + 58£ 3 Wi (16) + 1170£ 2 VFi (15) + 9660^1 (14) + 27300Wi (13) ) 
+Am (e 2 ^i (17) - ll^i {16) - 409^!^ 



+4o3 (s^W^ + llie^i (16) + 969 Wi< 15 >) . 

= A 2W + ^(16) + 504|l£ + 32090l£ W U + 188m3i ^(13) 



+2018016 Wi {12) 



(A6) 



+ (3 £4)30 + l/2£ 3 4n - 2 £4)i 2 ) Wi 17 + (54430 - 364i2 + f £ 2 4n) W4 {16) 
+ 8i3 W^ 15 ^A 1U + 2016 Wi( 14 )A m . 

= \A 2W (25 Wi (17) £ 5 + 1839 Wi (16) £ 4 + 50451 Wi (15) £ 3 + 641802 Wi (14) £ 2 + 3762486 Wi (13) £ 
+8072064 Wi (12) ) 
+i Am (W 17 ^ 3 + 51 W^t 2 + 813 W^ 15 ^ + 4032 W^ 14 ^ 
+3430 (Wi( 17 )£ + 18^i( 16 )) -2 4)12 (W 17 )£ + 18Wi( 16 )) . 



(A7) 



(A8) 
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= \A 2W (ihW^i 7 + 1523 W^ 16 \ G + 61761 WV 15 >£ 5 + 1289820 W^ 1A \ A + 14889420 Wi( 13 )£ 3 
+94316040 Wi {l2) i 2 + 300900600 Wi (11) £ + 369008640 W ni!; ' 



+2 Am (Vi (17) £ 5 + 76 Wi (16) £ 4 + 2156 Wi (15) £ 3 + 28392 Wi (14) £ 2 + 172536 W^C 

V (A9) 
+384384 Wi (12) J 

-A 012 (Wi^ 3 + 173 Wi^ 2 + 1475 VFi( 15 )e + 1024 VFi^ 14 )) 

3^030 (Wi (17) £ 3 + 1 15 Wi (16) £ 2 + 1249 Wi (15) £ + 3392 Wi (14) ) . 

= -|^3oo (W 17) £ 12 + 177 Wi (16) £ n + 13413 Wi (15) £ 10 + 572670 WV 14 )£ 9 + 15257970 W^ 13 ^ 8 
+265552560 W^ 12 ^ 7 + 3072429360 Wi (11) £ 6 + 23597814240 W^ 10 ^ 5 
+118118800800 Wi (9) £ 4 + 370767196800 W x (8) £ 3 + 681080400000 Wi {7) £ 2 
+642939897600 W x (6) £ + 228843014400 W\ (5) ) 
-3^201 (5 Wi (17) £ 10 + 709 W^i 9 + 42051 Wi^f 8 + 1365672 W^ u ^ 7 
+26708136 VtV 13) £ 6 + 325909584 Wi (12) £ 5 + 2487204720 W x (11) £ 4 
+11568997440 Wi (10) £ 3 + 30849698880 Wi (9) £ 2 
^41565363840 W x {i) i + 20704844160 Wi 



+^ 



/ ' (A10) 

-2 A120 ( 9 Wi (17) £ 8 + 976 Wi {m) i 7 + 42728 Wi (15) £ 6 + 978432 Wi (14) f 5 + 12689040 W^ 13 '^ 

+94174080 Wi (12) C 3 + 383423040 W^ 11 ^ 2 + 761080320 Wi (10) £ + 536215680 Wi (9) ) 

+2A102 (W 17) £ 8 + 108 W^ 16 ^ 7 + 4704 W"/ 15 ^ 6 + 107016 W^ 14 ^ 5 + 1375920 W^ 13 ^ 4 

+10090080 Wi (12) £ 3 + 40360320 W x (11) £ 2 + 77837760 W x (10) £ + 51891840 W x (9) ) 

-4m (5 Wi (17) ^ 6 + 393 Wi (16) C 5 + 11511 WiC (15) + 157140 Wi (14) C 3 + 1016820 Wi (13) £ 2 

+2784600 Wi (12) ^ + 2325960 Wi (11) ) 

+9A o3 (Vi (17) £ 6 + 77 Wi (16) ^ 5 + 2223 Wi (15) ^ 4 + 30180 Wi (14) £ 3 + 196980 Wi (13) £ 2 

+556920Wi (12) ^ + 491400 Wi (11) ) . 
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= -f A 300 (^ 10 W^ + 149£ 9 W^ 16 ) + 9339£ 8 W^ 15 ) + 322728£W 14 ) 

+6772584^ 6 VFi (13) + 89618256£ 5 VFi (12) + 751710960£ 4 VFi (11) + 3912068160£ 3 VFi (10) 

+11961069120£ 2 W^( 9 ) + 19148088960^1 ^ + 11987015040 W x (7) ) 
-A 201 (^W^ + ^ + + 329721£ 5 W 14 ) + 4672395^ ^i (13) 

+38640420£ 3 VFi( 12 ) + 180360180^ + 429188760^1 ( 10 ) + 392432040 

-A 120 (8^W^ 17) + 657C 5 VFi (16) + 20769£ 4 V^i (15) + 320964£ 3 ^i (14) + 2532348£ 2 VFi (13) (All) 

+9546264^1 ( 12 ) + 13189176VFi( n )) 
+3^102 (VWi (17) + 83^ 5 VFi( 16 ) + 2653£ 4 W^ 15 ) + 41468£ 3 tt^ 14 ) 

+330876C 2 VFi (13) + 1260168fWi (12) + 1753752VFi (11) ) 
-3^021 (? 4 VFi (17) + 71^ 3 V^i (16) + 1583£ 2 VFi (15) + 12858^i (14) + 29022VFi (13) ) 
+9^003 (^ A Wi^ + 57£ 3 VI^ 16 ) + 1093£ 2 W^ 15 ) + 8070^i( 14 ) + 17850^i( 13 )) . 
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o = ^a 300 (e 14 ^i (17) + 205£ 13 ^i( 16 ) + i82i5^ 12 )vFi( 15 ) + 925260c 11 w x ( 14) 

+29849820£ 10 VFi (13) + 642762120£ 9 Wi {l2) + 9454044600£ 8 W / i (11) + 95610715200£ 7 Wi 10 

+660712852800£ 6 WV 9 ) + 3062396937600C 5 ^i (8) + 9210023222400£ 4 

+17036091072000e 3 VFi (6) + 17762576832000£ 2 

+8935774848000^1 (4) + 1525620096000^1 (3) ) 
-3^201 (W 17 )<£ 12 + 169 W^ 16 ^ 11 + 12 157 W^ 15 ^ 10 + 489230 W^ 14 ^ 9 + 12178530 W^ 13 ^ 8 

+195839280 W^ 12 ^ 7 + 2063421360 Wi (11) £ 6 + 14153499360 W 1 (10) £ 5 

+61556695200 Wi (9) £ 4 + 161124163200 W x (8) £ 3 + 230659228800 W x (7) £ 2 

+148929580800 W 1 ( 6 >$ + 25427001600 Wi ( 5 )) 
-| A120 (2 W// 17 ^ 10 + 269 VFi^ 16 ^ 9 + 15009 VF^ 15 ^ 8 + 453768 W^ 14 \ 7 + 8148504 V^( 13 ^ 6 12) 

+89618256 W x (12) £ 5 + 600359760 Wi (11) £ 4 + 2355312960 Wi (10) £ 3 + 4955670720 W x (9) £ 2 

+4618373760 Wi (8) £ + 1089728640 W x (7) ) 
+ i^i 02 (w/ 17 ^ 10 + 133 W^ 16 ^ 9 + 7323 Wi^ 8 + 217896 W^ 14 ^ 7 + 3837288 W^ 13 ^ 6 

+41185872 Wi {l2) i 5 + 267387120 Wi (11) £ 4 + 1006125120 W x (10) £ 3 + 1997835840 Wi (9) £ 2 

+1712430720 W x (8) £ + 363242880 Wi {7) ) 
-A 02 i (2 + 191 W^ 16 ^ 7 + 7171 V^i^ 15 ^ 6 + 135846 W/ 14 ^ 5 + 1385370 Wi^ 4 

+7502040 wV 12 )f 3 + 19819800 W^ 11 ^ 2 + 20900880 W^ 10 ^ + 5045040 W^ 9 )) 
+3 A 003 (W 17) £ 8 + 97 W^ 16 ^ 7 + 3713 Wi^ 6 + 72090 W^ 14 ^ 5 + 759150 Wi^ 4 

+4291560 Wi (12) £ 3 + 12022920 Wi (11) £ 2 + 13693680 Wi (10) £ + 3603600 Wi (9) ) . 

= 1 ,4 2 i (10 Wi (17) £ 9 + 1297 Wi (16) £ 8 + 69655 Wi (15) £ 7 + 2022230 Wi (14) £ 6 + 34761090 Wi (13) £ 5 

+364236600 Wi (12) £ 4 + 2306424120 Wi (11) £ 3 + 8421613200 Wi (10) £ 2 

+15881065200 VFi (9) ^ + 11589177600 VFi (8) ) 
+3 Am (W 17) £ 7 + 101 VFi (16) £ 6 + 4067 VFi (15) ^ 5 + 84140 Wi (14) £ 4 + 959140 Wi (13) £ 3 

+5973240 Wi (12) £ 2 + 18618600 W x + 22102080 Wi (10) ) (A13) 
+9 A 030 (W 17 ^ 5 + 59 W^ 16 ^ 4 + 1227 W^ 15 ^ 3 + 10946 W^ 14 ^ 2 + 40838 W/ 13 ^ 

+53872 VFi (12) ) 

-3A i2 (w/ 17 ^ 5 +41 Wi^k 4 + 273 W^ 15 ^ 3 - 6146 W^ 14 ^ 2 - 78638 W^ 13 ^ 
-206752 Wi (12)N ) . 
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= \A 2 iq (jWx^H 11 + 1107 W^ 16 ^ 10 + 74133 W 1 ( 15 ^ 9 + 2760408 W x (14) £ 8 

+63115416 Wi^i 7 + 923792688 Wi (12) ^ 6 + 8757180432 Wi (11) £ 5 + 53215081920 (10) £ 4 
+200587907520 ( 9 ^ 3 + 439575776640 W x ( 8 ^ 2 + 493647073920 V^i (7) £ 
+209227898880 T^i (6) ) 
+2A U1 (W 17) £ 9 + 126 Wi (16) £ 8 + 6546 Wi (15) C 7 + 182868 Tyi (14) £ 6 + 3004092 I^i (13) £ 5 

+29811600 W^ 12 k 4 + 176576400 W^ 11 ^ 3 + 592431840 + 998917920 T4V 9) £ (A14) 

+622702080 VFi (8) ) 

+3^030 (W 17) £ 7 + 81 W! {16 k 6 + 2535Wi (15) C 5 + 39444 VFi (14) £ 4 + 329364 T^i (13) C 3 

+1500408 iyi (12) ^ 2 + 3583944 W^ ll) i + 3459456 WV 10) ) 
+A i2 (W 17) £ 7 + 129 WV 16) £ 6 + 6027 VFi (15) £ 5 + 134484 Wi (14) £ 4 + 1544004 W^ 13 ^ 3 

+8969688 iyi (12) C 2 + 23633064 W x + 20756736 W^A . 



Appendix B: Complex potentials and superintegrable systems in £7(1, 1) 
admitting separation in parabolic coordinates 

As a byproduct of this study, we have obtained a complex potential in E2 that can be 
viewed as real potential in the pseudo- Euclidean place E(l, 1). This is the only complex 
potential in Ref. 22] which is separable in parabolic coordinates and does not admit a 
Killing vector. We also give the subcase of this potential which admits a Killing vector and 
its third-order integrals. 



1. System admitting a single third-order integral:^ 

The following potential 



2 , , JX 2 I 'Ol - if ) , 7 



C + IT] i + IT) 

(3(2x + iy) 7 

= 2ax - 



2ar cos (6) 



\Jx + iy ^/x~~fTy 

r (2 cos 6 + i sin 6 1 ) /3 + 7 



a/ r (cos 9 + i sin 1 
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admits one third-order integral with A 3 = — A 02i , A i 2 = 2iA 02 i and the remained Ajke = 0. 
The functions G\ and G 2 are given by 

d = i (-2iria£ 2 + 4£ar] 2 - 2ir)£ p + 2iar) 3 + r) 2 fi + 7) A 021 
G 2 = (-2£aif + Aii 1 a£ 2 + 2irj^p + 2^a + f3 £ 2 -7) A 021 . 

Putting iy = t, we obtain a superintegrable system in E(l, 1); namely, 

2 y/X + t y/X + t 

2. System admitting group symmetry 

A special case of the previous potential admits a Killing vector. The potential is: 

Vl = j^—- (B3) 

The third-order integral is defined by the following non-zero Aj ke , along with an arbitrary 
constant k\ which is the coefficient of the Killing vector 

^4o2i ~\~ 2iA§\ 2 

Ax)3 — ^030 = 7, 1 ^120 — —iAiu + A102, Aiq 2 , Am, A021, ^012, 



Gi = 7^012 + il A 021 + i 7 (f + ir?) 2 A102 - I (C 2 - -4m + &i (f + "7) 
G 2 = i 7 A)i2 + 7 Am + I (f + A102 - y - ^) An + ik i (f + ^) • 
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